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Abstract 

A number of physical systems exhibit a particular form of asymptotic conformal invariance: 
within a particular domain of distances, they are characterized by a long-range conformal interac- 
tion (inverse square potential), the apparent absence of dimensional scales, and an SO(2,l) sym- 
metry algebra. Examples from molecular physics to black holes are provided and discussed within 
a unified treatment. When such systems are physically realized in the appropriate strong-coupling 
regime, the occurrence of quantum symmetry breaking is possible. This anomaly is revealed by 
the failure of the symmetry generators to close the algebra in a manner shown to be independent 
of the renormalization procedure. 
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I. INTRODUCTION 



An anomaly is defined as the symmetry breaking of a classical invariance at the quantum 
level. This intriguing phenomenon has played a crucial role in theoretical physics since its 
discovery in the 1960s [lj. In addition to its use in particle phenomenology of the standard 
model I2J and its extensions, it has been a fruitful tool for the study of conformal invariance 

n 

in string theory |3J. 

Surprisingly, the presence of an infinite number of degrees of freedom does not appear 
to be a prerequisite for the emergence of anomalies. This fact was first recognized within 
a model with conformal invariance: the two-dimensional contact interaction in quantum 
mechanics ^|. In conformal quantum mechanics, a physical system is classically invariant 
under the most general combination of the following time reparametrizations: time trans- 
lations, generated by the Hamiltonian H; scale transformations, generated by the dilation 
operator D = tH — (p • r + r ■ p) /4; and translations of reciprocal time, generated by the 
special conformal operator K = 2tD—t 2 H+mr 2 /2. These generators define the noncompact 
SO(2,l) « SL(2,R) Lie algebra U 

[D, H] = -ihH , [K, H] = -2ihD , [D, K] = ihK . (1) 

This symmetry algebra has also been recognized in the free nonrelativistic particle lal, the 
inverse square potential j?l Isj], the magnetic monopole jjj, the magnetic vortex jlo| . and 
various nonrelativistic quantum field theories 0, 111! Q] . Furthermore, conformal quantum 
mechanics has been fertile ground for the study of singular potentials and renormalization, 



using Hamiltonian [13, 



14 



15| as well as path integral methods Most importantly, a va- 
riety of physical realizations of conformal quantum mechanics have been recently identified, 
as discussed in the next section. 

The main goals of this paper are (i) to illustrate the relevance of conformal quantum 
mechanics for several physical problems, from molecular physics to black holes, and (ii) to 
show the details of the breakdown of the commutator algebra (JIJ for the long-range conformal 
interaction. In Sec. Owe introduce a number of examples that can be regarded as physical 
realizations of conformal quantum mechanics. In Sec. 11111 we show that the origin of the 
anomaly can be traced to the short-distance singular behavior of the conformal interaction. 
In Sec. |lV|we introduce a generic class of real-space regulators, within the philosophy of the 
effective-field theory program. In Sec.|Swe compute the anomaly for the regularized theory 



and show that it is independent of the details of the ultraviolet physics, and in Sec. IVII 
we comment on various renormalization frameworks. After the conclusions in Sec. I VI II we 
summarize a number of technical results: a derivation of the anisotropic generalization of 
the conformal long-range interaction (Appendix EJ ; a study of interdimensional dependence 
(Appendix [BJ; a derivation of the near-horizon properties of black holes (Appendix [CJ) ; and 
a derivation of useful integral identities ( Appendix ID]) . 

II. RELEVANT PHYSICAL APPLICATIONS 

In recent years, diverse examples of systems have been studied from the viewpoint of the 
conformal algebra (0), assumed to be a representation of an approximate symmetry within 
specific scale domains. In the applicable conformally invariant domain, the relevant physics 
is described by a d-dimensional effective Hamiltonian 



which involves a long-range conformal interaction; or, alternatively, by its anisotropic coun- 
terpart 

H=£--±F(fl), (3) 
2m r z 

where ft stands for the angular variables and F(fl) is a generic anisotropy factor that 
accounts for the angular dependence. Equation (j3J) is discussed in Appendix lAl 

In the problems considered below, A = 2mg/h 2 is the dimensionless form of the coupling 
constant and v — (d — 2)/2; furthermore, the choice h — 1 — 2m will be made for the 

n 

problem involving black holes. In all cases, the strong-coupling regime [1J] is defined by the 
condition g > g(*\ with a critical dimensionless coupling A**) = X\ — (I + v) 2 (for angular 
momentum I), when the Hamiltonian model (J2J) is adopted In addition, in the strong- 
coupling regime, as deduced in Sec. IIHI an uncontrolled oscillatory behavior of the Bessel 
functions of imaginary order iO makes the conformal system singular and regularization 
is called for. The characteristic parameter = ^/X — (I + v) 2 strictly corresponds to the 
Hamiltonian (J2J); in physical applications, such as those of Sees. Ill K\ 111 Bl and 111 Cl we will 
define 



H 



P 2 9 
2m r 2 ' 



(2) 




(4) 
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which will turn out to be crucial in parametrizing the anomalous physics of the conformal 
system in the presence of symmetry breaking. In discussing these realizations, we will 
explicitly use a subscript to emphasize the effective nature of the parameter of Eq. 
as arising from a reduction framework. The same notational convention will apply to the 
dimensionality (d c g). As shown in Appendix [BJ even when interdimensional equivalences 
are introduced, the value of the parameter (j3J) is a dimensional invariant. 

A. Dipole-bound anions and anisotropic conformal interaction 

The three-dimensional (d e g = 3 or h> c g = 1/2) interaction between an electron (charge 
Q = — e) and a polar molecule (dipole moment p) was the first physical application to 
be recognized as a realization of this anomaly [17[. When the molecule is modeled as a 
point dipole, this interaction can be effectively described with an anisotropic long-range 
conformal interaction of the form (JHJ): V(r) = —gcos8/r 2 , in which the polar angle 9 is 
subtended from the direction of the dipole moment. For this potential, the dimensionless 
coupling is A = — 2m K e pQ/h 2 = p/po, with m being the reduced mass of the system and 
K e the electrostatic constant. Thus, the relevant scale for phenomenological analyses is 
provided by po ~ 1.271 D (where D stands for the debye). 

As shown in Appendix El in some sense, the generic anisotropic conformal interac- 
tion © — of which the electron-molecule interaction is a particular case — can be reduced 
to an effective isotropic conformal interaction for the zero angular-momentum channel [see 
Eq. flA7|) ]: this corresponds to an effective Hamiltonian of the type (j2J), with an appropri- 
ate effective coupling A e ff. More precisely, this equivalence is achieved, after separation of 
variables in spherical coordinates, at the level of the radial equation. In addition, the cor- 
responding value of A e ff is identical to the eigenvalue 7 of the angular equation, which is a 
function of the dipole coupling A. The effective conformal parameter becomes 



where X *^ = u 2 s for each eigenvalue 7 of the ang ular equation. When this outline is im- 
plemented, according to the procedure of Ref. |l7( or its generalization of Appendix El the 
existence of a critical dipole moment p^ for binding is predicted; the order of magnitude of 
its "conformal value," A w « 1.279, or p^*> £s 1.625 .D, has been verified in numerous exper- 




(5) 
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iments jl^l [Q . In particular, when binding occurs, extended states known as dipole-bound 
anions are formed. These conclusions have also been confirmed by detailed ab initio calcula- 
~ Q and by studies that income the effect, of rotational de gre es of ffeedo, Q, 



tions 



18, 



which also modify slightly the value of p(*\ 

In short, the central issue in this analysis — also shared by the other physical realizations 
discussed in this paper — is the existence of a conformally invariant domain whose ultraviolet 
boundary leads to the anomalous emergence of bound states via renormalization. As a 
result, these states break the original conformal symmetry of the model and modify the 
commutators as we will show in the next few sections. This simple fact alone captures the 
essence of the observed critical dipole moment in polar molecules and leads to an analytical 
prediction for the energies of the conformal states, as discussed in Sec. I VII and Appendix IA1 



B. Near-horizon black hole physics 



A generic class of applications of conformal quantum mechanics arises from the near 



horizon conformal invariance of black holes, its impact on their thermodynamics 



2l|, and 



its extension to superconformal quantum mechanics [22]. In particu 



tlio Hamiltoiiian (2| have been used to explore horizon states 3, 24\ and to shed light on 
black hole thermodynamics j^J . Another class of current applications j3] involves a many- 
body generalization of Eq. (J2J): the Calogero model, which has also been directly linked to 
black holes 2(|. These remarkable connections seem to confirm the conjecture that it is the 
horizon itself that encodes the quantum properties of a black hole (27 1. 

In this context, we consider the spherically symmetric Reissner- Nordstrom geometry in 
D spacetime dimensions, whose metric 



ar, analyses based on 



ds 2 = -f(r) dt 2 + [/(r)]" 1 dr 2 + r 2 dtt D _ 2 



(6) 



is minimally coupled to a scalar field with action (c = 1 and h = 1) 



S 



-~ / d D x 



-g [g^ d M $ d v § + mV] 



(7) 



In Eq. ©, dQo-2 stands for the metric on the unit (D — 2)-sphere, f{r) = 1 — 2 (a_M-/r) ~ + 
(6 ( g/r) 2( - D ~ 3 ' ) , and the lengths a m and 6q are determined from the mass M and charge Q 
of the black hole respectively 28]. In this approach, the conformal structure is revealed 
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by a two-step procedure discussed in Appendix and consisting of: (a) a reduction to an 
effective Schrodinger-like equation, to be analyzed in its frequency (u) components; (b) the 
introduction of a near-horizon expansion in the variable x = r — r + [with r = r± being the 
roots of f(r) = 0, and r + > r_]. Two distinct scenarios emerge from this reduction: the 
extremal and nonextremal cases, when r + = r_ and r + 7^ r_, respectively. We will omit any 
discussion of the extremal case, which is known to pose a number of conceptual difficulties 
and is otherwise beyond the scope of the framework presented in this paper. As for the 
nonextremal case, the following facts arise from this reduction: 

(i) The ensuing effective problem is described by an interaction 

, . (near horizon) , , 

V(x) oc -x- 2 , (8) 

which is conformally invariant with respect to the near-horizon coordinate x. 

(ii) The effective Hamiltonian, still being a d- dimensional realization of the conformal 
interaction, does not have the usual form corresponding to the radial part of a multidimen- 
sional Schrodinger problem. In particular, the angular momentum variables appear at a 
higher order in the near-horizon expansion. 

(iii) The coupling constant A e fr is supercritical for all nonzero frequencies. This can be 
seen from Eq. ([Clip , which implies that 

The conclusion from this procedure is that the relevant physics occurs in the strong-coupling 
regime, in which the framework discussed in this paper can be applied. 

C. Other applications 

While Sees. Ill Al and 111 Bl conform to the title of this paper, applications in other areas 



29, 30] stands out. This effect 



of physics are also likely. Among these, the Efimov effect 
is expected to arise in a three-body system with short-range interactions, in which at least 
two of the two-body subsystems have virtual or bound s-states near zero energy. As in the 
case of the dipole-bound anions of Sec. Ill Al these are spatially extended and weakly bound 
states. Unfortunately, the combination of phenomenological parameters needed to form 
these states, together with their weakly bound nature, has defied experimental detection to 
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date. Nonetheless, this effect is regarded as relevant in the description of the three-body 
nucleon interaction [31]. The most outstanding feature of these three-body interactions in 
three spatial dimensions is the fact that these problems are reduced to an effective equation 
with a long-range conformal interaction in the strong-coupling regime. In terms of possible 
experimental detection, this effect is currently being studied for the description of various 
systems, including helium trimers and nuclear three-body halos [3^. 

The conformal nature of the effective interaction, for the three-body systems described 
above, can be deduced as follows. Typically, one starts by introducing hyperspherical coor- 
dinates with hyperradius p = r, in a c4ff-dimensional configuration space for the internal de- 
grees of freedom; if the one-particle dynamics occurs in a d- dimensional space, then ci e ff = 2d 
for the internal dynamics of the three-body system (as the total number of coordinates is 
3d, but d of them are eliminated in favor of the center-of-mass coordinates). Consequently, 
when a hyperspherical adiabatic expansion 3^| is combined with a Faddeev decomposition 
of the wave function a reduction to a d e & = 2d realization of our conformal model (J2J) 
is obtained. These conclusions can be gleaned from the conformal nature of the effective 
adiabatic potential VMr) arising ncn this .eduction f_t k Q, 

^ff(r) = -ff> A eff = (ci-l) 2 + 9 e 2 ff , d eS = 2d, (10) 

where g e ff and A e ff are related as described after Eq. (jSJ). 

Incidentally, due to the interdimensional equivalence of Appendix [Bl this result is often 
quoted in its one-dimensional reduced form [from Eq. (jB4|) ]. X(d = 1) = A c fj — (d— l) 2 + l/4 = 
2 ff + 1/4. For example, for the all-important case of ordinary three-dimensional space, 
d e ff = 6 and X e s = 4 + 2 ff . Furthermore, the coupling constant in Eq. (fTU|) depends 
upon the physical parameters defining the system: when the scattering lengths are large, 
it is function of the three ratios of particle masses. In particular, for the lowest angular 
eigenvalue of a three-body three-dimensional system of identical bosons with zero-range 



'ormal parameter (j3J) is approximately given 



two-particle interactions, the characteristic con 
by the solution of the transcendental equation 

8sinh(^=V3e eff cosh(^ , (11) 

so that O e ff ~ 1.006, which corresponds to the strong-coupling regime. 

In short, the essential feature shared by the problems discussed above is the existence 
of an effective description in terms of SO (2,1) conformal invariance, which results from a 



prescribed reduction framework. We now turn our attention to this generic effective problem, 
characterized by the Hamiltonian of Eq. (J2J). As different dimensionalities are required for 
the applications to which Eq. (J2J) refers, we will analyze this problem for the arbitrary d- 
dimensional case. Our goal is to investigate and characterize the possible realization of a 
conformal anomaly within this scope. 



III. CONFORMAL ANOMALY AND SHORT-DISTANCE PHYSICS 



Conformal symmetry is guaranteed at the quantum level when the naive scaling of op- 
erators, described by the algebra ((H), is maintained. A measure of the deviation from this 
scaling is afforded by the "anomaly" j3] 



.4(r) = ±[D,H] + H 



„d-2 



/'" " _ r V(r) 



r d-2 



V(r) (12) 

(13) 



2 

(valid for arbitrary d spatial dimensions), in which 1 is the identity operator and £ r = r • V. 
At first sight, the right-hand side of Eq. (fT2*|) appears to be zero for any scale- invariant 
potential; however, upon closer examination, this apparent cancellation may break down at 
r = 0, where the interaction is singular. Equations (|T2^) and (JTHJ) can be directly applied to 
any of the interactions within the conformal quantum mechanics class, and reduce to the 
familiar results known for the two-dimensional contact interaction j^, However, the 
most interesting case is provided by the Hamiltonian (j2J), whose symmetry breaking can be 
made apparent by means of the formal d- dimensional identity V- [r/r^ 1 ] = Qd-i S^(r), in 
which Qd-i is the surface area of the unit (d — l)-sphere S x ; then, 

^(r) = -s%V- 2 5W(r). (14) 

Despite its misleading appearance, this term is not identically equal to zero, due to the 
singular nature of the interaction at r = 0. The recognition of this remarkable singular 
term, as well as of its regularized and renormalized counterparts, leads to the central result 
of our paper: the proof of the existence of a conformal anomaly. 

However, two important points should be clarified. First, Eq. (|T4"j) is merely a formal 
identity, whose physical meaning can only be manifested through appropriate integral ex- 
pressions. Second, the coordinate singularity highlights the need to determine the behavior 
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of the wave function near r = 0. Therefore, nontrivial consequences of Eq. can only be 
displayed by the expectation value with a normalized state \^f), 

j t (D)„ = (A(r)) 9 = -9^y~ J d d r6^(r) |r^(r)| 2 . (15) 

A similar analysis applies to the anisotropic interaction of Eq. (J3J); in this case, 

j t (D) 9 = (A(t)) 9 = -9^ J d d r6^(r) |r^(r)| 2 F(O) . (16) 

It should be noticed that the intermediate steps leading to Eqs. (|T5Jl and (JTfij) are based on 
formal identities involving the (^-dimensional S function. For the unregularized inverse square 
potential, the integrals in Eqs. (|15j) and (|16j) select the limit r — » of the product r u ^f(r), 
which is known to be proportional to a Bessel function of order iO, with defined in Eq. Pj). 
This limit is ill defined in the strong-coupling regime, due to the uncontrolled oscillatory 
behavior of the Bessel functions of imaginary order. Consequently, a regularization procedure 
is called for; inter alia, this procedure will assign a meaningful value to Eqs. (fT5|) and (fTtjj). 



IV. REGULARIZATION AND RENORMALIZATION: THE EFFECTIVE-FIELD 
THEORY PROGRAM 

The Hamiltonian (j2J, in the strong-coupling regime, describes an effective system with 
singular behavior for short-distance scales. This interpretation, in which regularization 
and renormalization are mandatory, is inspired by the effective-field theory program 3(|. 
The required regularization procedure is implemented in real space, where the ultraviolet 
physics is replaced over length scales r < a. The effective theory that comes out of this 
renormalization is expected to be applicable within a domain of energies of magnitude \E\ <S 
E a = h 2 /2ma 2 . The scale E a defines an approximate limit of the conformal regime from 
the ultraviolet side; effectively, this limit prevents the singular interaction from yielding 
unphysical divergent results for supercritical coupling. 

Specifically, we consider a generic class of regularization schemes that explicitly modify 
the ultraviolet physics; each scheme is described by a potential V^^r), for r < a, where a 
is a small real-space regulator. An appropriate procedure for the selection of solutions of 
this singular conformal interaction was proposed in Ref. J^, using a constant potential for 
r < a. Our approach is based on a generalization of this method, in which a core interaction 
V^(r) is introduced. 



Incidentally, in this section, we consider a core V«)(r) = V«\r) with central symme- 
try V^(r). Even though this condition is not strictly necessary, it leads to a tractable 
derivation. Moreover, it is also consistent with the original rotational invariance of the 
isotropic singular interaction and captures the singular behavior of the potential, which 
originates from its radial dependence (even in the anisotropic case). The generalization for 
an anisotropic conformal interaction is nontrivial, but when this interaction is reduced to 
an effective radial problem, the procedure developed in this section can be applied. 

The core interaction is subject to the conditions of finiteness 

-oc<V = min[FM(r)]=-|^ (17) 

and continuous matching with the external inverse square potential at r = a, 

00(a) = V {> \a) = -g/a 2 . (18) 

It should be noticed that these restrictions imply that Vq < or K > 0, and that K = A + q, 
where £ > is the dimensionless energy difference between the minimum Vo and the matching 
value (|18|). In addition, in this approach, the energies for the interior problem will be 
conveniently redefined from the minimum value Vo; specifically, 

U(r) = V {< \r) -V , E = E-V . (19) 

For the spherically symmetric long-range conformal interaction of Eq. central sym- 
metry leads to the separable solution 

tt (r) = y «°»("M0 } [ dn^lYUm 2 = , (20) 



I — I 

in which Y/ m (0) stands for the ultraspherical harmonics on S^ -1 [38], which have been 
conveniently redefined with a normalization integral equal to the solid angle Qd-i- Then, 
the corresponding effective radial Schrodinger equation for v(r) becomes 



d 1 1 d 
dr 2 r dr 



2 



k - ^~t^ ~ V(r) 



v{r) = , (21) 



where V(r) = 2mV(r)/h 2 and k 2 = 2mE/h 2 . In particular, for bound states, k = in and 
Eq. (|21|) provides solutions of the form 

v (<)(r) = Bi jV wi +u (kr; k) for r < a , 

Mr) = S ( 22 ) 
v(>)(rj — A\ )V KiQ\Kr) for r > a , 
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in which Kie(z) is the Macdonald function [3J|, and where k is defined from E = h k : /2m, 



so that k = \/—k, 2 — Vq, with Vo = 2mVo/h 2 < 0. The regularizing core is arbitrary and 
Wi +U (kr; k) is a particular real solution in that region, 



d 2 1 d 
dr 2 r dr 



?-M_ M(r) 



w i+u (kr; k) = 



(23) 



where W(r) = V(r) — Vo; as an example, wi +u (kr; k) is a Bessel function of order I + z/ when 
the potential V^(r) is a constant. 

The solution f!22l) can be completely determined by enforcing the following three addi- 
tional physical conditions: (a) continuity at r = a of the radial wave function; (b) continuity 
at r = a of the logarithmic derivative of the radial wave function; and (c) normalization 
of the wave function. In what follows, these conditions will be stated using the auxiliary 
parameters 

f = «a , i=ka, (24) 

which satisfy Eq. (JTHJ), i.e., 

e 2 + e = n • (25) 

Consequently, these conditions (a)-(c) become, respectively, 



B itU wi +u {£; k) = A t)V Kie(g) 



(26) 
(27) 



and [cf. Eq. flUJ)] 



/poo 
d d T |$(r)| 2 = y drrKr)| 2 = 1 , (28) 

where we have conveniently redefined the logarithmic derivatives from C\q (£) = 
^[lniQe(0]' w hh £^ = £<9/<9£, and similarly for in terms of wi +1/ (^;k). Ex- 

plicitly, Eq. (|28p takes the form 



in which the normalization constants can be chosen to be real, and where 



and 



f£ r _ "I 2 

Ji+u(^,k)= / dz;z w i+u (z;k) 
Jo 1 i 



(29) 



(30) 



(31) 
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Equations and then provide the values of the constants Ai >v and Biy, for example, 

k) 



B, 



tj l+v {h~k)+ 



KM 



2 A -V2 



(32) 



For reasons that will become clear in the next section, it is convenient to rewrite Eqs. (|3U|) 
and (|31J) in an alternative way, using the generalized Lommel integrals of Appendix [DJ 
First, the integral defined by Eq. (|3()jl . which applies to the external domain (r > a), can be 
expressed as 

JC l& (0 = l[K t e(0] 2 M%\0, (33) 

where 

M%\Z)=[£&\Z)] 2 + &-e. (34) 

Similarly, the integral defined by Eq. 1)31 J), which applies to the internal domain (r < a), 
takes the form 

1r - - i2 . .» _ _ 

(35) 



where 



and 



M\t ] Mk)= c\ii{i-k) + e-d+vf-eu&k) 



•(<). 



r _ 1 2 

Ui +U (^;k)= dzz wi +u (z;k) 
Jo L 



1 + -^ )£/(*;*) 



(36) 



(37) 



with U = U/E and £ z = zd/dz. The Lommel integral relation ()33|) appears to be simpler 
than Eq. ()35|) because of the absence of an extra core U (z; k) in the external domain. 

In addition, the continuity conditions of the potential, Eq. (j!8|) . and of the logarithmic 
derivatives, Eq. (|27j) . imply the equality of the "matching functions" an d Q36p. i.e., 



(38) 



As a corollary, a combined Lommel relation can be obtained by elimination of the matching 
functions from Eqs. (jHSJ) and (JSSJ), 

2 



k)-U l+v (l;k) 



fc) 



(39) 



Even though the implementation of a renormalization procedure is a necessary condition 
for the emergence of the conformal anomaly, the actual details of this procedure are not 
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explicitly required. It suffices to know that these details are to be consistently derived from 
Eqs. ()24)) -()39 )) . which permit the exact evaluation of all relevant expectation values, and by 
enforcing the finiteness of a particular bound state energy. 



V. COMPUTATION OF THE CONFORMAL ANOMALY 



The value of the anomalous part of the commutator [D, H] is given as the "anomaly" 
A(r) in Eq. (|12jl. In Sec. this quantity was computed for the unregularized inverse 
square potential in terms of the formal identity (J 14)) ; this expression, in turn, led to an ill- 
defined expectation value ()15|) . This difficulty can be overcome when the singular conformal 
interaction is regularized according to the generic scheme introduced in Sec. IIV1 Then, 
Eq. ()12j) will in principle yield two different contributions: one for r < a and one for r > a, 
with the latter being of the form ()14j) : thus, 



AJr) 



1 + 2* 



6{a 



g 



r d - 2 5 {d \r)9{r 



(40) 



where 8(z) stands for the Heaviside function, is the regularized counterpart of Eq. (|14j). 
Explicitly, this leads to an expectation value 



d 



«) 



v(>) 



(41) 



where the integration range is split into the two regions: < r < a and r > a. Moreover, 
the identically vanishing second term 



(42) 



in Eq. (|4(J|) shows that the source of the conformal anomaly is confined to an arbitrarily small 
region about the origin. This result can be confirmed from a straightforward replacement 
of Eq. (JI2J) by A{r) = -\(d - 2) V(r) + \ V-{r V(r)}, which is identically equal to zero for 
any domain that excludes the origin, when applied to any homogeneous potential of degree 
-2 (a defining characteristic of the external conformal interaction). 

Once Eq. (|4*2"|) is established, the anomaly can be computed from the contribution arising 
from the ultraviolet domain r < a, 



| (D)* = (A(r))if = [ d d 

Jr<a 



l + -£ r )V{r 



(43) 
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In Eq. (J33J), V(r) = V (< )(r) can be replaced using Eq. (jTUjl . and \l/ a (r) using Eqs. (|50jl 
and (j22j) ; when these substitutions are made and the dimensionless variable £ in Eq. (|24|) is 
introduced, Eq. (|4~5j) becomes 



P 



c?2; z 



wi +v (z; k) 



1 2 



Vn + (7 - 



Despite its cumbersome appearance, the integral in Eq. 
the definitions (|3Tj) and (|37)l are introduced, so that 



(44) 



can be easily evaluated once 



E dt E k 2 



V Ji +v (£; k)+EU l+v (i; k) 



hJi+vfck)+ J l+V (£;k)-U l+U fck) 



(45) 
(46) 



where Vq was replaced through the relation (jl9j) or ()25|). and i£ = — /i 2 K 2 /2m. Furthermore, 
in Eq. (jjBj) , the difference Ji +U (^k) — Ui +V (£; k) can be evaluated employing Eq. so 
that 



--(D) 
E dt { h 



~ ~ -I 2 

fc) 



(47) 



Finally, the coefficient B^ u can be eliminated using Eq. (J32J), which shows that the right- 
hand side of Eq. (}4*Tj) is identically equal to one for any bound state. This remarkable 
simplification concludes the proof that the anomaly defined in Eq. (|12|) is indeed given by 

d 



dt 



(48) 



where E is the energy of the corresponding stationary normalized state. 

In short, we have validated the relation (J48)) — which agrees with the formal prediction 
from properties of expectation values 3^|. This validation has been established using a 
generic regularization procedure. Therefore, regardless of the renormalization framework 
used, an anomaly is generated. The generality of Eq. (}4~Sj) makes it available for a variety 
of physical applications, and is a necessary condition when the theory is renormalized. 



VI. RENORMALIZATION FRAMEWORKS 



In the previous section we showed that the property (J48)) and related symmetry-breaking 
results are independent of the details of the regularization procedure. Because of the gener- 
ality of the real-space regularization approach presented in this paper, these results extend 



the two-dimensional analysis of Ref. 



40j in a number of nontrivial ways: 
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(i) For arbitrary renormalization frameworks, other than the "intrinsic" one of Ref. 
(see below). 

n 

(ii) For any dimensionality d. Again, the two-dimensional case of Ref. 401 has unique 
features that considerably simplify the derivation within the intrinsic framework. This is 
particularly relevant because the physical applications that appear to be most interesting 
are (^-dimensional realizations of this phenomenon, with d = d e g 7^ 2. 

(iii) For any bound state and angular momentum channel (and not just for the I = 
channel associated with the ground state considered in Refs. 34, 35, 4oT|). 

In this section we highlight the relevance of these results with an overview of the real-space 
"effective," "intrinsic," and "core" renormalization frameworks (according to the presenta- 
tion of Ref. j^). and discuss their relationship to the present anomaly calculation. Despite 
their apparent differences, these frameworks share the basic physical requirement that the 
system is renormalized under the assumption that the ultraviolet physics dictates the pos- 
sible existence of bound states of finite energy; the corresponding energies E and values of 
k oc y[^[ are then required to remain finite. 

In order to facilitate the implementation of this renormalization program, it is convenient 
to display the specific limiting form that Eq. (|27jl takes when a — ► 0; more precisely, 

cot [a, (6, Ka)\ ^ , (49) 

where C«\ka;k) (k ~ <1} £«>(N) and 

a (9, no) = 6 [in ) + 7e ] , (50) 

with 70 = — {phase [r(l + z©)]} /© (which reduces to the Euler-Mascheroni constant 7 
in the limit — > 0). 

In the effective renormalization framework, the system is regularized maintaining finite 
values of \E\ E a == h 2 /2ma 2 . This condition defines an asymptotic conformally invariant 
domain; within that domain, the condition Ka < 1 limits the ultraviolet applicability of 
this effective scheme. Most importantly, this condition is systematically applied to derive 
physical predictions in a direct manner, within the prescriptions of Sec. IIVI As a result, 
Eq. (fUJj) leads to the bound-state energy levels J^lJ 

E n = E exp ( — — ) , (51) 
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in which Eq < is an arbitrary proportionality constant. This derivation also shows that, 
as ultraviolet physics sets in for \E\ > E a (that is, for Ka > 1), no claim can be made as to 
the nature of the states on these deeper scales. 

A few comments are in order regarding Eq. (|51|) . First, it explicitly displays a breakdown 
of the conformal symmetry, by the introduction of a scale \E \ and an associated sequence 
of bound states. Second, the scale \Eq\ arises from the renormalization procedure. Third, 
as a renormalization scale, |-Eb| cannot be predicted by the conformal model and it is to 
be adjusted experimentally. Fourth, once the experimental determination is carried out, an 
unambiguous prediction [from Eqs. (jIHJ) and (|5Tj) ] follows, 

E n+ i j t / 27r\ 

" ex P -7T > ( 52 ) 



dt 

within the range of applicability, sa C 1. This is in agreement with the conclusions of 
phenomenological analyses of the Efimov effect JscJ. 

The alternative intrinsic and core frameworks are characterized by the fact that the limit 
£ = Ka — > is strictly applied before drawing any conclusions about the physics. Therefore, 
in order to keep the bound-state energies and k values finite, a running coupling parameter 
is explicitly introduced, so that Eq. ()49|) is still maintained in this limit. The running 
parameter is either the conformal coupling g, in the intrinsic framework, or the strength N 
of the regularizing core interaction, in the core framework. 

In the case of the intrinsic framework, the dependence g = g(a), equivalent to = 9(a), 
is enforced. This leads to the asymptotic running behavior ~ 0, which ensures that the 
left-hand side of Eq. (J49|) remains well defined. This limiting procedure leads to the renormal- 



ization framework of Refs. 



13, 



14 1 ; in particular, Eq. (|52|) [with the condition ~ 0] implies 



the existence of a single bound state. In its original form, the renormalization framework of 



Refs. 







14 1 was based upon a Dirichlet boundary condition, which we now reinterpret as 



an effective Dirichlet boundary condition 16] u(r = a) ^ 0, for the reduced radial wave 
function u(r) = y/rv(r). This result is guaranteed by the prefactor y/r, regardless of the 
behavior of v(r). As for v(r), two distinct cases should be considered: (i) the special case 
characterized by the simultaneous assignments d = 2, I = 0, and constant V^(t) [or, to be 
more precise, with q = |Vo|a 2 — A = o(a 2 )], for which £ = ka = 0(0) and cosa(0, Ka) ^~°' ) 0; 
(ii) the generic case, characterized by d ^ 2, or I ^ 0, or V^(r) not being constant, for 
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which the variable £ = ka acquires a nonvanishing limit value = y/Jl + + i [as either 
(/ + z/) 7^ or £ 7^ 0], and sina(6, Ka) 0. 

The smallness of the variable £ for case (i) above is the main reason for the simplicity of the 
derivation of Ref. j^j]. In effect, in this case, Eq. (J43)) can be approximated using the small- 
argument behavior of Bessel functions without explicitly computing full-fledged Lommel 
integrals. Thus, Ji +U=0 (ka) (a ~ 0) 6 2 /2 and \Bi =0)U=0 \ = \Ai =QtU=0 Kie(na) / J (ka)\ (a ~ 0) 
«/(V7re), leading to d jdt (a ~ 0) 2ttBJ JqM^oA 2 (a ~° } E, as discussed in Ref. Q- 
By contrast, for the generic case (ii), the analysis presented in this paper, based on the 
theory developed in Sec. IIVI and Appendix lD| is inescapable. 

Finally, in the core renormalization framework, the strength of the core interaction be- 
comes a run ning coupling parameter: N = N(a), but the conformal coupling g remains 



constant As a result, Eq. (jlUj) provides the limit 



in renormalization analyses of the three-body problem 



eye 
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31 



e running that has been used 
Ac 



Incidentally, the "effective" renormalization framework discussed in Ref. |4l( (and sum- 
marized in this section) leads directly to a characterization of the thermodynamics of black 
holes. In essence, this amounts to a reinterpretation of 't Hooft's brick wall method 22(, in 
which ultraviolet "new" physics sets in within a distance of the order of the Planck scale 
from the horizon. The computation of Appendix IU1 shows that the leading behavior near the 
horizon is conformal and nontrivial, in that the effective system is placed in the supercritical 
regime. This asymptotic leading contribution, governed by the effective conformal inter- 
action, requires renormalization and provides the correct thermodynamics [44]. It should 
be noticed that there is an alternative treatment, based upon the method of self-adjoint 



extensions, which has been recently discussed in Refs. (23 



24|. 



VII. CONCLUSIONS 



Realizations of the conformal anomaly involve a breakdown of the associated SO(2,l) 
algebra. In this paper we have shown that the actual emergence and value of the conformal 
anomaly rely upon the application of a renormalization procedure, but are otherwise inde- 
pendent of the details of the ultraviolet physics. In this sense, the results derived herein are 
robust and totally general. As such, they are intended to shed light on the physics of any 
system with a conformally invariant domain for which the short-distance physics dictates 
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the existence of bound states. 

In particular, the dipole-bound anions of molecular physics and the Efimov effect are phys- 
ical realizations of this unusual anomaly. In addition, the intriguing near-horizon physics 
of black holes appears to suggest yet another example of this ubiquitous phenomenon; the 
details of the thermodynamics arising from this conformal description will be reported else- 
where. 
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APPENDIX A: ANISOTROPIC LONG-RANGE CONFORMAL INTERACTION 
AND CONFORMAL BEHAVIOR OF DIPOLE-BOUND ANIONS 

In this appendix we show the mathematical procedure that reduces the anisotropic inverse 
square potential to an effective isotropic interaction. 

The Schrodinger equation for the Hamiltonian (jSJ) can be separated in spherical coordi- 
nates by means of 



As a result, the angular part S(O) of the wave function is no longer a solution to Laplace's 
equation on the unit (d — l)-sphere S d_1 ; instead, it satisfies the modified equation 




(Al) 



with normalization 




(A2) 




(A3) 



where 



A = -A 2 + XF(Q) 



(A4) 
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and A 2 = L 2 /h 2 is the dimensionless squared angular momentum. The corresponding radial 
equation 

' A; + u(r) = (A5) 



dr 2 \ r 2 

is coupled to Eq. (|A3J) through the separation constant 7. Equation (|A5|) can be compared 
against the radial equation of an isotropic inverse square potential, which is obtained by 
another Liouville transformation j3] of the form (jAlj) . but with ultraspherical harmonics 
instead of H(O) and for Kff(r) oc r~ 2 without angular dependence; the effective equation 



d 2 u(r) 



+ 



k 2 + 



A eff - (I + v? + 1/4 



u{r) = (A6) 



<ir 2 

is identical to Eq. (|A5|) when the following identifications are made: 

Kff(r) = -^f , 9eS = ^^eff , A eff | /=0 = 7 ■ ( A ?) 

Consequently, Eq. (jA5|) can be thought of as the radial part of a rf-dimensional effective 
isotropic conformal interaction for 1 = 0. 

Furthermore, the values 7 are quantized from the angular equation (jA3|) and depend upon 
the coupling A of the anisotropic potential. This relationship can be made more explicit by 
expanding, in the ultraspherical-harmonic basis Yi m (Q), the anisotropy factor 

F(n) = J2 F i™Yi m (n) (as) 

Z,m 

and the angular wave function 

H(ft) = ]TH /m >WO). (A9) 

This decomposition yields the matrix counterpart of Eq. (jA3|) . whence the anticipated rela- 
tionship can be formally displayed by the infinite secular determinant 

£>(7,A) = detM(7,A), M( 7 , A) = -A(X) + 7 1 , (A10) 

in which 1 is the identity matrix; the matrix elements in Eq. (|A10|) are 

(/m|M( 7 , A)|ZW> = [1(1 + 2u) + 7] 6 u ,6 mm , -\J2 Tl 

m,Z'm';Z"m" ^"m" , (AH) 

l",m" 

where 

Wm';/»m» = J dfirf-l Y^ft^'m" {Sl)Y Vm , (fi) . (A12) 
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Finally, the components S/ m of the angular wave function can be formally obtained for every 
eigenvalue 7 in the usual way, and satisfy [from Eq. (|A2|) ] 



= |2 



1 . 



(A13) 



Lm 



As an example of this general theory, one can consider the particular three-dimensional 
case (v = 1/2) of the electron-polar molecule interaction described in Sec. Ill Al In this case, 
the matrix elements (|A11|) become 



(lm\M(j,X)\l'm!) = [1(1 + 1) + 7] 6 U , 5 mm > — X\\j ^™^L ™\ ftv-i 



(l + m + l)(l-m + l) 

+ \ I —■, . TTTTT, . —\ <V 7+1 mm t 



(A14) 



(21 + l)(2/ + 3) 

which correspond to a matrix of block-diagonal form with respect to m and tridiagonal in 
I. Then, the secular determinant (jAlOj) factors out in the form D("y,\) = n m D m (7,A), 
with the reduced determinant -D m (7, A) in the m sector; thus, for given m, the equation 
det M (7, A) = implies that 



£>m(7,A) 



7 



(2 + 7) 







'15 



(6 + 7) 



. 



(A15) 



Equation (jA15|) has been used for the determination of the critical dipole moment A^*-* ~ 
1.279 Q| when 7 = ■y^ . When the determinant is expanded (to high orders), additional 
roots appear for the critical condition 7^) = 1/4 and for different values of m. This pattern 
also illustrates how one would completely solve the generic anisotropic problem: Eq. (|A15|) 
or its generalization (jA10|) can be used to obtain the eigenvalues 7 that correspond to a 
given coupling A; these eigenvalues replace the usual angular momentum numbers. In the 
molecular physics case described above, the values of 7 can be easily evaluated numerically. 
When A < \(*\ no such values produce binding; a first "binding eigenvalue" 70.0 is obtained 
when A > \^ = 1/4, for the first root with m = 0; as the strength A of the interaction 
increases, a second binding eigenvalue 70,1 is produced for the first root with m — 1, when 
A ~ 7.58 or p rs 9.63 D; the next eigenvalue 7^0 arises from the second root with m = 0; 
etc. For each one of these values of 7 = 7j, m , an energy spectrum of conformal states is 
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governed by Eg . (I5T1), with e fj given by Eq. (J3j). These bound states have been observed 
experimentally isi for the case when 70,0 is the only binding eigenvalue, a condition 
that corresponds to typical molecular dipole moments. 

Most importantly, this analysis confirms that the conformal anisotropic problem can be 
reduced to the isotropic one, and the same symmetry-breaking considerations apply. 

APPENDIX B: DIMENSIONALITIES AND INTERDIMENSIONAL DEPEN- 
DENCE 

The spatial dimensionality d e g of a physical realization of conformal quantum mechanics 
is best characterized or defined as the dimension of the configuration space needed for 
a complete description of the dynamics within the conformal approximation. Typically, 
this quantity can be directly identified from the nature of the radial variable used in the 
description of scale and conformal symmetries. 

For instance, with this convention, molecular anions can be naturally seen as a three- 
dimensional realization (d e g = 3); the Efimov effect, in a <i-dimensional one-particle space, 
as a (2<i)-dimensional realization (<i e fr = 2<i); and the near- horizon conformal physics of black 
holes, in D = d + 1 spacetime dimensions, as a (^-dimensional realization (d e s = d). 

Of course, there is a certain degree of arbitrariness in the selection of d c g, due to the 
existence of a formal relationship connecting problems of different dimensionalities. This 
can be seen from the reduced Schrodinger-like radial equation of a conformal problem (J2J), 



d 2 u(r) 
dr 2 



, A - (I + vf + 1/4' 



k' + 



u{r) = . (Bl) 



Equation (jBl|) depends on the number of spatial dimensions only through the combination 
I + v, a property known as interdimensional dependence . As a consequence, the radial 
part of the solutions for any two conformal problems are identical when their coupling 
constants are related by 

X(d'; I') = \{d; /) + (/'- 1 + d'jlf -(1-1 + djlf . (B2) 

Moreover, 

e(d') = Q{d) (B3) 

is a dimensional invariant of these formal transformations. Correspondingly, the conformal 
physics is totally determined by the invariant value of this parameter. 
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However, the interdimensional equivalence of Eq. (jB2|) is severely limited by the fact 
that the full-fledged solutions (wave functions) are not identical, because the angular mo- 
menta are different in different dimensionalities. The only exception to this is the formal 
equivalence among the I = angular momentum channels of problems with arbitrary di- 
mensionalities (as these channels do not involve additional dimension-dependent angular 
variables); in particular, an effective one- dimensional coupling can always be introduced for 
a d- dimensional problem with I = 0: 

X(d' = = 0) = X(d;l = 0) + \- . (B4 ) 

Even in the special case of the equivalence described by Eq. (|B4[) . the full-fledged wave 
functions still retain a trace of the "physical dimensionality" d, because (with an obvious 
notation) u(r) = \E'| (J=1 (r) = r^ _1 ^' 2l i r |d(r); for example, in the case of the three-dimensional 
Efimov effect, the full-fledged wave functions are of the form \l/(r) oc r~ 5//2 w(r), reflecting 
the fact that d e{i = 6. 

The example of the near-horizon conformal behavior of black holes presents a number of 
peculiar features that deserve a separate treatment in Appendix 



APPENDIX C: NEAR-HORIZON CONFORMAL BEHAVIOR OF BLACK 
HOLES 

In this appendix we present an algebraic derivation of the conformal invariance exhibited 
near the horizon of a black hole. 

From Eqs. (JBJ) and (J7J), it follows that the equation of motion satisfied by the scalar field 
in the black-hole gravitational background is 

(□ - m 2 ) $ = -?=<9m (-yf^g cT d u <5>) - m 2 $ 

= -i$ + + (V + (jP ~ 2)/ ) $' + ^A D _ 2 $ - m 2 $ = , (CI) 

where the dots stand for time derivatives and the primes for radial derivatives in the chosen 
coordinate description of the background, while is the Laplacian on the unit (D — 2)- 

sphere. In addition, by separation of the time and angular variables, 

r, n) = e-^ip lm (r)Y lm (n) , (C2) 
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Eq. (|C1|) turns into 



u 2 m 2 



a 



^ r) + l '£ + C£^y (r)+ ^ 



<^(r) = 



(C3) 



with a = 1(1 + -D — 3) being the eigenvalue of the operator —A D _ 2 . Equation (jC3|) can be 
further reduced, by means of a Liouville transformation j3| 



<p(r) = g(r) u(r), g(r) = exp 
to its normal or canonical form 



/ 7" 



1/2 -(D-2)/2 



w"(r) + J(r) u(r) 



with normal invariant 

T/ . uj 2 m 2 



(D-2)(D-4) a 
4 / 



1 _ ir if _ p-2)r 



2 / 4/ 2 



2r/ 



(C4) 



(C5) 



(C6) 



The conformal behavior of the Schrodinger-like equation ()C5|) near the horizon can be 
studied by means of an expansion in the variable 



x = r — r i 



(C7) 



with r = r + being the largest root of f(r) = 0. The nonextremal case is characterized by 
the condition 

f+ = f'(r+) ± , (C8) 

equivalent to r + ^ r_. Then, 



f(r) =f' + x[l + 0(x)\ , f(r) =/;[! + 0(x)] , /"(r) = /J [1 + 0(*)] 



(C9) 



where /" = f"(r + ). Thus, with corrective multiplicative factors of the order [1 + 0(x)}, 
it follows that /"// ~ fl/(fLx) and /'// ~ while r ~ r + , so that the only leading 
terms in Eq. (Jd are u 2 / f 2 ~ cj 2 /(/» 2 and /' 2 /(4/ 2 ) ~ l/(4x 2 ). As a result, Eq. (Jd is 
asymptotically reduced to the conformally invariant form 



u"(x) 



(/+ 



x" 2 [1 + 0(x)]u(a;) = 



(CIO) 



where, by abuse of notation, we have replaced u(r) by Equation (|C10|) indicates the 

existence of an asymptotic conformal symmetry driven by the effective interaction 



A 



cff 



X- 



Aefr = ^ + e 2 ff , e 2 



cff 



L/'(r+)J 



(Cll) 
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as follows by rewriting Eq. (|C10|) in the d-dimensional format of Eq. (|B1|) . This proves the 
claims made in Sec. Ill Bl and, in particular, Eqs. (jHJ) and (jHJ). 

A final remark is in order. The effective Hamiltonian (jC10|) did not fall "automatically" 
within the d- dimensional format of Eq. (jBl|) . The extra terms —[(I + v) 2 — l/4)]/r 2 , usually 
obtained by reduction of a multidimensional Schrodinger equation in flat space, are still 
present, but at higher orders in the expansion with respect to the near-horizon coordinate 
x; in Eq. (jC6|) . they correspond to 



(£>-2)(£>-4) a 
4 / 



1 



(* + ") 
/ 



4 + ^ 



-7 



i 



-,=o\- 



x 



(C12) 



[with v = (d — 2)/2 = (D — 3)/2]. Thus, the angular momentum-together with its associated 
dimensionality variable — decouples from the conformal interaction (jCll|) in the near-horizon 
limit. It should be noticed that we had to rewrite Eq. (jC10|) in the 1 = 0, <i-dimensional 
format in order to present this problem within our unified conformal model (0). Alterna- 
tively, one could write Eq. (|C11|) in a simpler one-dimensional reduced form [from Eq. (jB4J) ]. 
X(d = 1) = A c ff — v 2 + 1/4 = Q 2 S + 1/4, with the same value for the dimensional invariant 

e eff . 



APPENDIX D: GENERALIZED LOMMEL INTEGRALS 



In this appendix we derive a generalization of the Lommel integrals [47[ for an arbitrary 
Sturm-Liouville problem 



L x v(x) = fig(x)v(x) 



d_ 

dx 



dx 



+ q(x) 



(Dl) 
(D2) 



and apply it to the reduced radial Schrodinger equation (|2Tj) . These generalized integrals 
are needed for the exact evaluation of expectation values in the anomaly calculation. 
In what follows, we rewrite the differential equation (jDl|) in the form 

^ \p{x)v'{x)\ = - [a 2 g(x) + q(x)] v(x) , (D3) 

with an eigenvalue fi = a 2 and where the prime stands for a derivative with respect to x; 
moreover, v (x) can be chosen to be a real function. Next, after conveniently multiplying 
both sides by 2p(x)v' (x) , and integrating them with respect to x, Eq. (jD3|) turns into 



[p{x)v'{x)\ 



■t'l> 



J-2 



*1 



dxp(x) \a 2 g{x) + q(x)~\ — [v(x)]" 

(ajJu 
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(D4) 



in which both the lower (x\) and upper limits (£2) are completely arbitrary. Finally, after 
integration by parts and rearrangement of terms, Eq. (jD4|) leads to 



a 2 / dx [p(x)g(x)]' [v (x)} = [v(x)Y 



p(x 



v'(x) 



v\x) 



+ p(x) [a 2 g(x) + q(x)] 



■'■■2 



j- 2 



dx \p(x)q(x)}' [v(x)Y 



(D5) 



which generalizes the well-known second Lommel integral |47| of the theory of Bessel func- 
tions. A similar procedure could be applied for a generalization of the first Lommel integral, 
but this is not needed for the present purposes. 

The integral relation (jD5|) can be rewritten in a convenient form for the reduced radial 
Schrodinger equation (|21|). which is of the generalized Bessel form 



{ dx 2 x dx ^ ^ ^ } 



v (x) = 



(D6) 



This is a particular case of the Sturm-Liouville equation (jD3|) . with density function g(x) = 
x, p(x) = x, and q(x) = — xW(x); however, it is also true that a straightforward set of two 
Liouville transformations ^| makes Eqs. (jDl|l and ()P6j) formally equivalent to each other. 
For Eq. ()P6|) . [p(x)q(x)]' = — [a; 2 W(x)]', and the final term in Eq. (|P5j) can be evaluated 
with the help of 

L [W{x) x 2 ] = 2x + ~ E x ) W{x) , (D7) 

where 1 is the identity operator and £ x = xd/dx, as in Sec. |lVl As a consequence, Eq. f)P5|) 
becomes 



f 2 1 

a 2 / dxx[v{x)} 2 = - [v{x)] 2 {{C{x)} 2 + [{ax) 2 - x 2 W{x)] } 



X-2 



XI 



XI 



(l.r.r[r(x)) 2 (t + ^E x ) >V(.r) 



;ps) 



where C(x) = xv'(x)/v(x) and both limits are still arbitrary. Equation (jP8|) is the desired 
generalization that can be directly applied to the reduced Schrodinger equations (J2l|) and 
to derive Eqs. and (|3K|). as we will show next. 
First, for the interior problem (r < a), Eq. ()P8j) turns into Eq. (jSSJ), by means of the 
substitutions 



x = r, a = k, x 2 W(x) = (I + v) 2 + r 2 lA{r) , v[x) = Wi +U (kr; k) , z = kr , (D9) 
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and with integration interval z G [0, £], where C, = ka. For this case, when r 2 U(r) — > 0, that 
is, for regular core potentials, the behavior of the differential equation at the origin implies 
that the contribution from the first term on the right-hand side of Eq. (jD8|) is zero for r = 0. 

Second, in a similar manner, for the exterior problem (r > a), Eq. (jD8|) turns into 
Eq. (jSHJ), by means of the substitutions 

x — r , a = k = in , x 2 W(x) = (I + u) 2 — A = — O 2 , v(x) = Ki®(nr) , z = nr , 

(D10) 

and with integration interval z G [£, oo], with ^ = «a. Here, the behavior of the differential 
equation at infinity implies that the contribution from the first term on the right-hand side 
of Eq. (jD8|) is also zero at that point. 
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